Quenching Factor for Low Energy Nuclear Recoils in a Plastic Scintillator 
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Plastic scintillators are widely used in industry, medicine and scientific research, including nuclear 
and particle physics. Although one of their most common applications is in neutron detection, 
experimental data on their response to low-energy nuclear recoils are scarce. Here, the relative 
scintillation efficiency for neutron-induced nuclear recoils in a polystyrene-based plastic scintillator 
(UPS-923A) is presented, exploring recoil energies between 125 keV and 850 keV. Monte Carlo 
simulations, incorporating light collection efficiency and energy resolution effects, are used to gen- 
erate neutron scattering spectra which are matched to observed distributions of scintillation signals 
to parameterise the energy- dependent quenching factor. At energies above 300 keV the depen- 
dence is reasonably described using the semi-empirical formulation of Birks and a kB factor of 
(0.014±0.002) gMeV~^cm~^ has been determined. Below that energy the measured quenching 
factor falls more steeply than predicted by the Birks formalism. 

PACS numbers: 32.50.-(-d; 78.70.Ps; 29.40.Mc; 28.20.Cz 



I. INTRODUCTION 

The response of organic scintillators to particle interac- 
tions in terms of the dependence on material, incident 
particle type and incident particle energy were first dis- 
cussed by Birks fH, . In general, it is found that the 
response arising from nuclear recoils (such as when irra- 
diated by neutrons) is significantly diminished in com- 
parison to the light output obtained from electron recoils 
(such as when irradiated by 7-rays). At higher energies 
(MeV and above), the scintillation output is generally 
found to be proportional to energy deposition but, at 
lower energies, a strong departure from proportionality 
has been observed for nuclear recoils. A thorough char- 
acterisation and understanding of such effects is essen- 
tial for accurate low energy calibration, especially given 
the widespread use of scintillators in contemporary sci- 
ence. One specific example, where the low energy re- 
sponse to nuclear recoils is paramount, can be found in 
the field of direct dark matter searches, both for the re- 
sponse of the dark matter targets themselves (e.g. noble 
liquid scintillators), and where scintillators find their ap- 
phcation in anti-coincidence detector systems It 
is in this context that the present results have been ob- 
tained. Conceptual designs for future, large active neu- 
tron rejection systems featuring scintillators are under 
discussion Q, and will require improved knowledge of 
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the low energy response, even when the main neutron 
detection mechanism is via radiative capture. In the case 
of polystyrene-based scintillators, little data exist for re- 
coils below ~1 MeV, which are produced, for example, 
by radioactivity neutrons. 

The present measurements were performed with 
the plastic scintillator used in the veto detector of 
the ZEPLIN-III dark matter experiment, based at 
the Boulby Mine, UK. ZEPLIN-III is a two-phase 
(gas/liquid) xenon detector designed to observe low- 
energy nuclear recoils from galactic weakly interacting 
massive particles (WIMPs) j8|-fl3|. For its second sci- 
ence run the detector has been enclosed by a polystyrene, 
(CsIIg)^, based veto detector The veto instrument 
includes 52 modules individually coupled to photomul- 
tiplicr tubes (PMTs), totalling ~1 tonne of scintillator, 
which, in the form of a barrel and a roof, surround the 
WIMP target. For a detailed discussion of the realised 
performance of the veto detector, the reader is referred 
to Ref. [13,]. 



A. Quenching 

The scintillation light yield for a nuclear recoil of a given 
energy is quenched, i.e., reduced in comparison to the 
scintillation output observed from an electron recoil of 
the same energy. A significant contribution to this differ- 
ence may be identified with the heat associated with the 
atom cascades generated by nuclear recoils as described 
by Lindhard |14| . A formalism in which the scintillation 
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light yield of highly ionising particles depends not only 
on the energy of the particle but also on its stopping 
power in specific materials was developed by Birks ll | , 
of which a detailed description is presented in Ref. |15[. 
The resulting relation may be written as: 
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where dL/dr is the scintillation yield per unit path length 
r, iS" is the absolute scintillation factor, BdE/dr is the 
density of excitation centres along the recoil ionisation 
track and A: is a quenching factor. By finding the ratio be- 
tween the light yield for electron recoils, Le, and for ions, 
ii, Eqn. ([l} may be rewritten in terms of the quenching 
factor for nuclear recoils, Qi, in integrated form as: 
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From ([U and ([2]), an energy dependence of the quench- 
ing factor is apparent. This is especially significant for 
the low energy region where the stopping power experi- 
ences greatest variation. The majority of the measure- 
ments obtained to-date for the quenching factor in plastic 
scintillators concentrate on neutrons and protons in the 
energy region above ^^1 MeV ;16i-i20i] . In recent years, the 
need for precise knowledge of neutron quenching factors 
for materials used in the direct search for dark matter 
has led to significant new measurements at low energies, 
often making use of dedicated neutron scattering facili- 
ties [2l],[22|. However, no recent measurements have been 
reported for plastic scintillators despite their incorpora- 
tion into several low energy event experiments. In this 
paper we present measurements of nuclear recoil quench- 
ing factors for energies below 1 MeV down to 125 keV. 



II. EXPERIMENTAL SET-UP 

One of the 52 plastic scintillator modules of the 
ZEPLIN-III veto detector was used for data taking in 
the Boulby Underground Laboratory, an intrinsically low 
background environment. The scintillator bar has the 
form of a parallelepiped of length 1 m, width 15 cm, and 
a trapezoidal cross-section with parallel sides of length 
15 cm and 12 cm. The polystyrene-based scintillator 
material (UPS-923A, p-terphenyl 2%, POPOP 0.02%, 
produced by Amcrys-H, Kharkov, Ukraine i23j]) has a 
density of 1.06 g/cm^ and a refractive index of 1.52 for 
blue light. The average light output for electron recoils 
has been measured to be ~5500 photons/MeV [3]. The 
scintillation light shows a peak intensity at 420 nm, a rise 
time of 0.9 ns and a decay time of 3.3 ns. The average 
bulk attenuation length for the 52 modules has been ex- 
perimentally measured and is found to be approximately 

1 m g. 



To increase the effective attenuation length of the plas- 
tic and improve light collection, a specularly reflective 
aluminised Mylar foil is placed at one end. Additionally, 
the module has been wrapped in diffuse reflector PTFE 
sheet on all sides. Light produced in the scintillator is 
detected with a PMT (ETEL-9302KB) of quantum effi- 
ciency 30% fsll optically coupled to the end opposite the 
mirror. 

Energy spectra were recorded with the dedicated data 
acquisition system of the veto detector (CAEN model 
V1724), which digitises waveforms with 14-bit resolution, 
0-2.25 V input range, 40 MHz bandwidth and a sam- 
pling rate of 10 MS/s. In this instance waveforms were 
parameterised using a bespoke data reduction software 
adapted from that developed for the ZEPLIN-III instru- 
ment [23 | . The trigger was provided by an external pulse 
generator at a constant frequency. Additionally, during 
the neutron source measurements, data were taken si- 
multaneously with a single-channel pulse height analyser 
("MAESTRO SCA"), triggered by an internal discrimi- 
nator. 

To measure the response to nuclear recoils, the scin- 
tillator was exposed to neutrons from a ^'*^Am-Be (a,n) 
source and, separately, to a ^^^Cf fission source. The 
plastic was shielded from 7-ray emission from the sources 
and the environment by enclosing it in a 20 cm thick cas- 
tle composed of low-background Cu and Pb in equal parts 
with an additional 4 cm of lead on the roof. Neutron ex- 
posures were performed with the sources placed directly 
on the castle (~50 cm above the sealed scintillator). Sys- 
tematic uncertainties in the setup were explored exten- 
sively from which it was found that variation in neutron 
source position had negligible effect. 

Crucially, 7-ray attenuation and external electron- 
recoil contamination within the nuclear recoil data have 
been quantified using Monte Carlo simulations and ded- 
icated measurements (see Sec. HI for detailed discussion 
of simulations). In particular, we examined the effect of 
varying the thickness of lead shielding placed over the 
castle. 7-ray emission spectra from the ^^^Cf and Am-Be 
sources (reconstructed from values given in the NuDat 
database [25}) have been studied separately. The ac- 
tual 7-ray activities were 21,000±2,100 7/s for the ^^scf 
source and 6,300±400 7/s for the Am-Be source (the lat- 
ter accounts only for the two highest energy 7-rays of 
3.21 MeV and 4.44 MeV from de-excitation of ^^C* pop- 
ulated by the Be(Q;,n) reaction). The simulations indi- 
cate that a single 7-ray from the Am-Be source would 
be transmitted through the shielding along with every 
30,000 neutrons (of which ^600 deposited energy in the 
scintillator bar) for the nominal lead thickness in our con- 
figuration, while no 7-rays from the ^^^Cf source expo- 
sure would be observed. Thus, the results show that the 
7-ray fluxes from the sources make no significant contri- 
bution to the neutron exposure data. To confirm this 
conclusion, an extended exposure of the scintillator to a 
11 kBq ™Co 7-ray source (1.17 and 1.33 MeV 7-rays) 
placed externally on the upper surface of the enclosure 
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was performed. No measurable increase in event rate over 
background was observed. Given that contributions from 
the 7-rays coming from the sources themselves are negli- 
gible in the neutron measurements, most 7-rays detected 
during the neutron exposure are generated internally (in- 
elastic scattering and radiative neutron capture). Non 
source-related backgrounds, arising, for example, from 
low level activity of shielding components and the plas- 
tic scintillator itself, are measurable but not significant 
above a threshold of 2 photoelectrons (phe) . 



III. SIMULATIONS 



The methodology used to extract the quenching fac- 
tor was first applied to liquid argon scintillation by the 
WARP group [2^; other examples followed d, 
Experimental data are compared to a comprehensive 
Monte Carlo simulation that includes a detailed descrip- 
tion of the experiment. The relationship between real 
energy deposition and resulting scintillation production, 
i.e., the energy-dependent quenching factor, is included 
as a parameter in the simulation. An iterative process 
is used to optimise the quenching factor, minimising on 
in the comparison between data and simulated energy 
spectra. The simulations have been performed with the 
GEANT4 toolkit ^Qj (version 9.2, with neutron cross- 
sections from ENDF/B-VI [3l|) using standard neutron 
spectra for the two sources (Am-Be ISO 8529-1 [13, ^^^Cf 
fission spectrum from S0URCES-4C [33,]). Emitted neu- 
trons and their secondaries are propagated including all 
relevant nuclear and electromagnetic physical processes; 
a set of optical processes describes the generation and 
detection of scintillation light from nuclear and electron 
recoil interactions in the scintillator. These photons are 
tracked to the photocathode of the PMT including rel- 
evant optical effects (reflection, refraction, attenuation) 
at which point the production of photoelectrons is sim- 
ulated. Appropriate random fluctuations are included 
to model the production of scintillation photons and the 
production of photoelectrons from the PMT photocath- 
ode. 

It should be noted that for a full description of neu- 
tron scattering in hydrogenous materials the standard 
GEANT4 elastic scattering process must be supple- 
mented with a model (G4NeutronIIPThermalScattering) 
to describe the energy region below 4 eV for the correct 
treatment of thermal neutron scattering from chemically- 
bound atoms. In these molecules several temperature- 
dependent vibrational modes are possible, which alter 
the scattering cross-section @| . This is of particular rel- 
evance to this study and radiative capture on hydrogen 
is enhanced by '^20% over the standard model. 



IV. 7-RAY CALIBRATIONS 

By deflnition, the response of the plastic scintillator to 
7-rays is unquenched, allowing standard 7-ray sources 
to be used to determine the overall gain of the sys- 
tem. Moreover, it is expected that the GEANT4 simula- 
tions should provide an excellent match to the 7-ray cal- 
ibration data, validating most processes included in the 
physics model and the accuracy of the geometry imple- 
mentation. The PMT gain was set such that both single 
photoelectron (SPE) peaks and Compton edge features 
could be resolved in all spectra, allowing presentation of 
the data in terms of absolute numbers of photoelectrons. 

With the roof of the shielding castle open, calibration 
measurements with a ^'^''Cs 7-ray source (4.7 kBq) were 
performed. Figure [T] shows the acquired spectra in com- 
parison to Monte Carlo simulations. Data were acquired 
with the CAEN acquisition system (solid black spectra) 
with a trigger provided by an external pulse generator 
operating at constant frequency. Signal pulses were then 
extracted from the recorded waveforms. The result of 
a GEANT4 simulation of this exposure is shown by the 
dashed red line; excellent agreement across the full en- 
ergy range is demonstrated. The scintillator module used 
has an attenuation length of 80 cm and the photoelec- 
tron yield with the calibration source above the centre 
of the plastic (48.7 cm from the photocathode face) is 
measured to be ~44 phe/MeV. 

The present simulations do not include spurious ef- 
fects such as dark emission from the photocathode, after- 
pulsing, or 13^ radiation from '^^K contamination in the 
glass PMT envelope. Each of these effects are known to 
contribute at the single to few photoelectron level with 
signiflcant rate [l3|, [35|. Consequently, we impose a 5 
photoelectron analysis threshold on the 7-ray calibration 
data and, therefore, on our neutron scattering analysis 
and results. 
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FIG. 1: Energy spectrum acquired from a ^'^^Cs 7-ray source 
exposure. The data acquired with the CAEN acquisition 
(solid black spectrum) with a threshold of 5 photoelectrons 
is shown in comparison to the simulation data (red dashed 
spectrum). 
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V. NEUTRON EXPOSURES 
A. Nuclear recoils 

Data were accrued for a live time of 600 s from separate 
exposures to the Am-Be source (5,500±300 neutrons/s) 
and the ^^^Cf source (3,400±170 neutrons/s). Placing 
the sources externally to the copper-lead enclosure at- 
tenuates the 7-ray emission from the sources to a negli- 
gible level. The impact of the enclosure on the neutron 
fluxes is illustrated in Figs. [2] and [3] for the two sources. 
The figures show the neutron emission spectra, the en- 
ergy spectra as they enter the scintillator (both referring 
to the y-axis on the left), and the resulting nuclear re- 
coil energy depositions in the polystyrene (y-axis on the 
right). The spectra at the scintillator interface include 
single neutrons being recorded multiple times as they 
are scattered out of the scintillator and re-enter again 
after interacting with the shielding. The recorded en- 
ergy depositions are the total integrated signal from each 
individual neutron-induced recoil event. The shielding 
attenuates significantly the neutron flux, and scattering 
reduces the energies of surviving particles. Since this is 
a large effect, we quantified how the uncertainty in the 
lead thickness affects the neutron spectrum at the scin- 
tillator interface. Variations up to ±0.5 cm are found to 
be statistically insignificant. 

The impact of thresholds in the simulated neutron 
source spectra (50 keV in both instances) has been exam- 
ined. Reasonable extrapolations down to keV do not 
change the recoil spectrum above threshold much and the 
ensuing quenching factor analysis is affected very little. 
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FIG. 2: Energy depositions in the scintillator from neutron- 
induced nuclear recoils coming from an Am-Be source (red 
hatched spectrum - referring to the scale on the right). The 
y-axis on the left refers to the neutron flux from the source 
(black dashed spectrum) and the differential neutron spec- 
trum when entering the scintillator bar (black solid spec- 
trum). 
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FIG. 3: Energy depositions in the scintillator from neutron- 
induced nuclear recoils coming from an ^^^Cf source (red 
hatched spectrum - referring to the right hand scale). The 
left hand scale refers to the original neutron spectrum (black 
dashed spectrum) in comparison with the differential rate of 
neutrons entering the scintillator bar (black solid spectrum). 



B. Quenching factor 

Where data do not exist for quenching factors at low en- 
ergies, it is customary to assume an energy-independent 
quenching as determined at higher energies. Various con- 
stant quenching factors have been considered and then 
compared to the present experimental data. Figures |4] 
and[S]show the data from the Am-Be and ^^^Cf source ex- 
posures in comparison to simulations which assumed an 
energy independent quenching factor, with Qi—0.1 yield- 
ing, in both cases, the best fit to the measured data. For 
such a value, the nuclear recoil spectrum is quenched suf- 
ficiently such that the (un-quenched) peak at 2.218 McV 
from 7-ray emission following radiative capture of neu- 
trons on hydrogen can be resolved. This feature, appear- 
ing at ~90 phe (with a ^30 phe) may, thus, be used to 
normalise the energy scales of simulated to observed spec- 
tra and extract a quenching factor for the nuclear recoils. 
Both figures show that by adopting energy independent 
quenching factors, a discrepancy occurs below ~35 phe. 
Above this value the goodness-of-fit is determined by sta- 
tistical fiuctuations only in both cases. The data shown 
in Figs. Hand [5] were recorded with the MAESTRO SCA 
for the reason of better statistics at the position of the 
hydrogen capture peak. Subsequent analysis was mainly 
performed using data acquired with the CAEN system 
to avoid bias from threshold dependent trigger setups. 
Aside from counting statistics, the two recordings do not 
differ from each other at higher energies. 

At very low photoelectron values (^20) greater diver- 
gence is observed between the the Monte Carlo and the 
measured data (see insets in Figs. H] and [SJ indicating an 
energy dependent behaviour of the quenching factor at 
low recoil energies. The methodology used to derive this 
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FIG. 4: Background-corrected energy spectrum originating 
from irradiation with an Am- Be source (grey shaded area) in 
comparison with simulations using the quenching factor Qi 
as a constant parameter for the whole energy range. The 
best agreement with the real data is met by the curve fea- 
turing Qi^O.l (blue solid spectrum). The peak at ~90 phe 
is the 2.2 MeV radiative capture 7-ray from hydrogen. The 
inset shows the impact of different constant quenching fac- 
tors at low photoelectron values. A marked discrepancy be- 
tween simulation and data suggests that an energy-dependent 
quenching factor may provide a better physical description for 
low recoil energies. 



energy dependent behaviour is as follows: a hypotheti- 
cal Qi{E) function is composed from 14 values of recoil 
energy (125, 150, 175, 200, 225, 250, 300, 350, 400, 450, 
550, 650, 750, 850 keV) and interpolated linearly between 
these points; a constant behaviour is assumed below and 
above this range. Above 1 MeV, low statistics and the de- 
creasing gradient of the quenching factor preclude more 
in-depth analysis. For each combination of Qi{E) pa- 
rameters (from a limited grid, guided to cover reasonable 
ranges), the full simulation is performed and calcu- 
lated for the resulting match to the data. Below 5 phe, 
spontaneous SPE emission and other effects described in 
Sec. IV can make a significant contribution to the exper- 
imental data and, therefore, this region is excluded from 
the minimisation. The Qi{E) parameters are modified 
for each iteration until no significant improvement in 
can be obtained. 

Figure |6] shows the resulting energy- dependent quench- 
ing factor from minimising the overall foi' both 
datasets. Here the 5 phe analysis threshold allows mea- 
surements down to a nuclear recoil energy of approxi- 
mately 125 keV. In the sub-threshold region below 5 phe 
an even stronger decrease in the quenching factor with 
energy would be inferred from uncorrected data. The 
68% confidence intervals shown are determined by the 
envelope of regions built up from quenching factor model 
curves which fulfil the criterion of XmodeL< Xmin + Qy^ 
where Qy= 15.89 for 14 free parameters ^Sq. 

Figure [7] compares the ^^^Cf data (black hatched his- 



FIG. 5: Background-corrected energy spectrum originating 
from irradiation with a '^^■^Cf source (grey shaded area) in 
comparison with simulations using the quenching factor Qi 
as a constant parameter for the whole energy range and a 
close up of the very low energy part of the spectrum as an 
inset at the top right. 
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FIG. 6: The nuclear recoil quenching factor in polystyrene 
plastic scintillator UPS-923A as a function of recoil energy 
deposition extracted by mean of x^ minimisation from com- 
parison of simulations to data from a ^^^Cf (solid green) and 
an Am-Be (dashed blue) exposure, respectively. The hatched 
areas represent the 68% C.L. bands (^^^Cf /, Am-Be \ ). 



togram) with the best fit of the energy-dependent (red 
solid) and the best fit of the energy-independent simula- 
tion (blue dashed), respectively. The inset provides the 
same comparison but for the Am-Be study. 



C. Birks factor, kB 

Following the discussion in Ref. [T^, the absolute value 
of the quenching factor for specific materials is expected 
to depend only on the so-called 'Birks factor', kB^ inde- 
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FIG. 7: Simulations using an energy-dependent value for 
Qi{E) in comparison with background corrected data ac- 
quired with the CAEN ADC (black hatched spectrum) from 
irradiation of the scintillator with a ^^^Cf and Am-Be source 
(inset), respectively. The best fit using minimisation is 
shown by the red solid histogram (□). For comparison, the 
blue dot-dash spectrum (Q) shows the use of a constant 
quenching factor (from best fit to data). 



pendently of the particle type. Consequently, the relative 
scintillation yield curve may be estimated by incorporat- 
ing the appropriate energy-dependent stopping power for 
the specific particle species. The kB factor is then deter- 
mined by fitting Eqn. ([2]) to experimental data. 

At higher energies, contributions to the observed en- 
ergy depositions come predominately from the scattering 
of protons in the plastic scintillator. For lower energy de- 
positions, it is found that carbon nuclei (99% ""^^C) con- 
tribute over 30% of the overall nuclear recoil energy de- 
positions. This relative fraction rises almost linearly in 
the lower energy region reaching ~50% below 20 keV, as 
shown in Fig. [51 

The energy-dependent quenching factors derived here 
from the two neutron sources are in good agreement with 
each other and may therefore be combined. This is signif- 
icant, since the neutron spectrum from an Am-Be source 
is somewhat uncertain below a few hundred keV, al- 
though this is especially so for stronger sources than the 
one used here [37|. The ^^^Cf fission spectrum, which 
is known more precisely, yields very similar results. A 
combination of the two results, following the prescrip- 
tion for asymmetric errors in Ref. [ssl], is presented in 
Fig. [5] as the black solid line, with uncertainty repre- 
sented by the shaded band. The quenching factor is seen 
to have a significant energy dependence, increasing in 
gradient towards low energies. In general, the observed 
dependence is reasonably similar to that expected from 
the Birks formalism above about 300 keV, but it departs 
from the expected behaviour at lower energies. Fitting 
the present results in the range of 300 keV to 850 keV 
results in a kB factor of (0.014±0.002) gMeV-^cm^^ 
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FIG. 8: Fraction of nuclear recoil energy depositions coming 
from carbon nuclei relative to the proton recoil contributions 
in the plastic scintillator averaged from exposures to both 
Am-Be and ^^^Cf neutron sources. 



The error given is statistical only. This is also shown in 
Fig. [9l with the contributions from protons, from carbon 
ions and the sum shown separately. Stopping powers for 
protons and carbon have been taken from NIST [s^ and 
the SRIM Stopping Range Tables [i^], respectively. 
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FIG. 9: A fit of semi-empirical calculations of Birks for com- 
bined proton and carbon stopping powers from varying the 
kB factor (solid red) to the measured quenching factors (black 
with hatched error band) above 300 keV nuclear recoil energy 
yields : kB = 0.0135 gMeV~^cm~^ Additionally, curves as- 
suming scattering off protons or off carbon nuclei only are also 
shown. Below this energy a clear divergence of measurement 
from the Birks description can be observed. 



The kB factor resulting from fitting the present data 
to the Birks formalism above 300 keV may be compared 
with a previous value of kB=0.009 gMeV~^cm~^ re- 
ported for a-particle interactions in polystyrene-based 
plastic scintillator ([l5| and references therein). The 
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level of agreement is good, considering choice of data 
acquisition alone can produce discrepancies of a factor of 
two [IB] . The current results exhibit a slightly steeper de- 
pendence than expected from the Birks formalism, but 
interestingly, the same feature is apparent in all previ- 
ous measurement for organic scintillators presented in the 
above reference. Above about 300 keV, the present data 
broadly support the semi-empirical description of Birks. 

Below that nuclear recoil energy, a clear deviation from 
Birks is evident in Fig. ^ indicating that the fraction 
of scintillation generated by low energy nuclear recoils 
appears to decrease even more rapidly. As mentioned 
above, the analysis reported here has been limited to 
above 5 phe, to avoid complications that might be in- 
troduced by single photoelectron level processes not in- 
cluded in the simulations. However, not only would in- 
clusion of these effects increase the discrepancy further, 
but examination of the 3-5 phe region indicates the trend 
continuing, with an even stronger dependence. A physi- 
cal mechanism responsible for this behaviour is unclear. 
This is the first measurement to report on quenching fac- 
tors at these energies for polystyrene. 

VI. CONCLUSION 

Despite the common use of plastic scintillators in indus- 
trial and scientific applications, little experimental data 
exist for the correlation between nuclear recoil energy 
deposition and scintillation output, especially below en- 
ergies of 1 MeV. Consequently, where plastic scintillators 
are used in low energy applications, a constant quench- 
ing of nuclear recoils is often assumed. We have measured 
the energy-dependent quenching factor for nuclear recoils 
in a polystyrene-based plastic scintillator (UPS-923A) 
for recoil energies between 125 keV and 850 keV. The 
analysis is based on comparison of observations of nu- 
clear recoil spectra obtained with broad-band neutron 
sources with Monte Carlo simulations using the GEANT4 
toolkit. Critical to this methodology is the accuracy of 
the Monte Carlo simulations; these demonstrated excel- 



lent reproduction of a 7-ray calibration source down to 
the analysis threshold of 5 photoelectrons. Significantly, 
the energy-dependent quenching factor for nuclear recoils 
was determined from measurements made with two dif- 
ferent neutron source spectra, yielding the same result. 

We find that the Birks model describes reasonably the 
relation between energy deposition and non-radiative 
transfer processes over part of the energy range studied. 
A Birks factor A;B=(0.014±0.002) gMeV-^cm-^ was 
extracted from the best fit between semi-empirical calcu- 
lations from a combination of proton and carbon nuclear 
recoils and the quenching factor curves presented. At 
lower energies a significant discrepancy between the 
Birks model and the present results was observed. 
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